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I. INTRODUCTION 



Traditionally the physics of the interaction of a nucleon with a nucleus has been rep- 
resented by an optical potential. Once that potential is specified, and by its use also the 
scattering matrix, all observables may be calculated. For intermediate energies there have 
been a number of methods of calculating the optical potential, the most successful of which 
have been the global Dirac phenomenology (DP) [1] and a microscopic Schrodinger model 
in coordinate space identified in the literature as the g-folding method. 

The global DP approach is the outgrowth of the seminal work by Walecka and his group 
with a second order reduction of the Dirac equation leading to Schrodinger equivalent equa- 
tions. Potentials in those equations have been found with which good fits to elastic nucleon- 
nucleus scattering data result |l( and more recently, using a global DP, very good fits have 
been found for nucleon energies ranging from 20 MeV to 1040 MeV jij. Such potentials are 
utilitarian for use in analyses of other data. In the g-folding approach, optical potentials 
are obtained by first deriving an effective NN interaction in-medium from a realistic free 
nucleon-nucleon (NN) potential. That effective interaction is then folded with a suitable 
representation of the density matrix of the target nucleus. 

Both approaches have had success. Not only have they given good predictions of differ- 
ential cross sections but also they have reproduced spin and integral observables 0, 0, IH El ■ 
A specific example of the success with the g-folding method has been its use in analyz- 
ing 200 MeV data from 208 Pb. Excellent results were obtained with structure given by 
a Skyrme-Hartree-Fock (SHF) model [tJ for the ground state that is consistent with the 
Friedman-Pandaharipande neutron equation of state 0. The global DP approach, based 
on the relativistic impulse approximation (RIA), fits the data well when proton and neu- 
tron densities, their root-mean-square radii, and the neutron skin thickness, all are credible. 
Targets of 40 Ca, 48 Ca and 208 Pb were considered 0. 

However, the two models fundamentally are different. The global DP model is dependent 
on the fitting of data to determine the parameters in the assumed potentials. Conversely, 
the microscopic Schrodinger model is dependent upon obtaining realistic effective NN inter- 
actions and upon the quality of the model chosen for the structure of the ground state of the 
target. A complex non-local optical potential results. Of course one can find a local equiv- 
alent to this but there is no a priori reason to assume that the (non-relativistic) potential 
deduced from the Dirac reduction resembles, in any way, that equivalent local potential. 

Yet, to date, only one real comparison jlOj] has been made simultaneously of the results 
from the Dirac and the ^-folding approaches for nucleon projectiles, and that for the 208 Pb 
nucleus and 100, 200, and 300 MeV nucleons. Also, an earlier comparison of the Dirac 
and the then favored phenomenological Schrodinger models revealed problems inherent in 
the latter The phenomenological Schrodinger model usually assumes a simple local 
Woods-Saxon form of the optical potential. Such an approach often fails to reproduce 
spin observables well enough due to assumptions made in defining the spin-orbit potential. 
However, a more complex phenomenological (local) Schrodinger model has been developed 
recently 11], with a global parameter specification, that does replicate very well a large data 
set for energies to 200 MeV. 

The past successes of the two model prescriptions suggest that both may be preset in a 
way that allows a fair comparison between them, particularly when the results are predictions 
in that no a posteriori adjustments to any facet of the base model details are considered. 
With that in mind, the prime purpose of this paper is to compare in detail results obtained 
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from the global DP model with those from the (/-folding model for proton scattering. First 
we report on results of predictions of angular observables for elastic scattering of 65 MeV 
and of 200 MeV protons from 12 C, 16 0, 40 Ca, 90 Zr, and 208 Pb. For all cases, the results of 
the calculations made were plotted before any data was added. We also consider integral 
observables for both proton and neutron scattering for energies to 800 MeV to evaluate the 
isospin dependence of the potentials. 

In this paper, the formalism of the microscopic (/-folding and that of the global DP 
approach are discussed briefly in Sect. II and the results are presented in Sect. III. In the 
latter there are two subsections. In the first of those, the results of calculations at 65 MeV 
and 200 MeV of differential cross sections and of spin observables, A y and Q, are compared 
with data. In the second, predictions of reaction and total reaction cross sections made 
using both methods are compared with data for the five nuclei selected to span the mass 
range. Conclusions that may be drawn are given thereafter in Sect. IV. 



II. MODELS OF THE NUCLEON-NUCLEUS OPTICAL POTENTIAL 

As detailed presentations of the two model prescriptions of interest have been published 0, 
only salient features are given herein. 



A. The microscopic (/-folding model 

To form the (/-folding model optical potential it is assumed that pairwise interactions 
between the projectile nucleon and each bound nucleon as prescribed by a large basis model 
of the target suffice. Thus, in the (/-folding approach []|, no phenomenological form for the 
potential is assumed nor is that potential derived from any other phenomenological potential. 
Instead the method begins with the (/-matrices of a realistic NN potential; those (/-matrices 
being solutions of the Brueckner-Bethe-Goldstone equations in infinite nuclear matter, 

g l s J(p>,p;k,K,k F ) = V3 T (p, P ') + -£ rV^ T (p',Q) 

in which Q(q, K, kf) is an angle-averaged Pauli operator with an average center-of-mass 
momentum K. The energies, E, in the denominator are single particle values within an 
averaged field 0. Such equations are solved for 32 NN channels, for Fermi momenta {kp) 
up to densities 1.5 times the nuclear density, and for many energies (k found from E with 
relativistic kinematics) to 300 MeV using the Bonn-B NN potential [12] as the starting 
interaction for the calculations. 

To form optical potentials in coordinate space then necessitates a mapping of those g- 
matrices onto Bessel transforms of a complex NN effective interaction. That has been made 
with particular emphasis placed upon reproduction of the on-shell values of the (/-matrices. 
The effective interaction is the mix of sets of Yukawa function form factors for central, 
two-nucleon spin-orbit, and tensor operators that can be used in the code DWBA98 (l3| . 
The density of the input structure for the target ground state is used to select what values 
of the complex and energy dependent strengths of the effective interaction are used in the 
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folding process. When folded with an appropriate ground state density of the target, the 
microscopic optical potential is obtained naturally, incorporating Pauli blocking and density 
dependences. The (coordinate space) potential contains both direct and exchange parts; 
the latter arising from antisymmetrization of the projectile and bound state nucleon wave 
functions. Consequently the potential is fully nonlocal. As formulated 0, there are no 
parameters in the model to adjust a posteriori. All results then are obtained from a single 
(predictive) calculation. For energies at which this theory is applicable, much success has 
been achieved predicting the observables from proton-nucleus scattering for a number of 
nuclei 0. 

Formally, and ignoring complications of spin interactions and expectations, we seek solu- 
tions of the Schrodinger equations of the form 



^V 2 - V c (r) + E 



*(r) = J U{r,r';E)^{r') dr' 



(2) 



where r, r' are relative NA coordinates, V c (r) is a Coulomb interaction and U(r, r') is the 
optical potential. In practice, such potentials, or their underlying multipoles, are not evalu- 
ated specifically when DWBA98 is used. However, it is useful to consider U(r, r') and how 
it may be formed to note how the target structure enters the process. Consider the overlaps 
which lead to U(r, r'), i.e. 

A 



U pA {0,l)= U pa (l,2,...,A) 



n=l 



<M1,2,...,A) =(*„|i^ e// (01)|* 



(3) 

where '0' denotes the projectile coordinates, and as all nucleons in the target are equivalent, 
we have chosen a specific entry ('1'). By so doing we can make use of a cof actor expansion 
of the nuclear many-body (ground) state, 



tf(l,2,... ,A)) = -Ly^ am (l)) q am |fr(l,2, 



A)) . 



(4) 



Here a specifies the quantum number set {n, l,j, (} where ( is the isospin projection. Then, 
as the underlined factor in Eq. (J3J) is independent of coordinate '1', 



|*) (V ? a'm'(l)|fl'e//(10){|^a m (l)) - | (0))} 



(5) 



when the required antisymmetry with projectile and struck nucleon is taken into account. 

The many-body reduced matrix elements of the particle-hole operator pairs are the OB- 
DME jij. Usually for elastic scattering they are just the nucleon shell occupancies of the 
target, r\ a . Thus such occupancies, the single nucleon bound (SP) states, and the effective 
interactions g e //(01), define the optical potential in this approach, as 



U{r x ,r 2 ;E) = ^Va 



(s) U D (R ls ,E) ip am (s) ds 



+ ^ am {v x )U Ex {R 12 ,E) Voim {v 2 ) 



(6) 



where R\ 2 = \r% — r 2 \, and U D and U Ex are appropriate combinations of components of 
the effective interaction for the direct and exchange contributions to the optical potential 
respectively 0. 
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B. Global DP optical model potentials 



A second-order reduction of the Dirac equation for nucleon-nucleus leads to a Schrodinger- 
equivalent equation that has physically correct, effective, central and spin-orbit potentials. 
Symmetry allows one to have Lorentz scalar, S, Lorentz vector, V, and tensor, T poten- 
tials which, collectively, are termed SVT potentials. However, the tensor force can always 
be replaced in a phase equivalent way by an effective relativistic potential of just scalar 
and vector parts, and as our approach is phenomenological, we consider the global DP us- 
ing scalar-vector (SV) potentials, together with the Coulomb potential. Solutions of those 
Schrodinger-equivalent equations have given S'-matrices with which observables have been 
predicted very well In the SV, the real vector potential is large and repulsive while the 
real scalar potential is somewhat larger and attractive. The imaginary vector potential is 
attractive and the imaginary scalar is repulsive. The Dirac equation so structured is suitable 
for simultaneous analyses of proton-nucleus and neutron-nucleus scattering data for incident 
energies up to several GeV. 

The scalar-vector (SV) model of global DP that has been fit to the elastic proton-nucleus 
scattering observables Q takes the form 

U(r,E,A) = V v (E,A)r(r,E,A) + V S (E, A) f s (r, E, A) 

+ iW v (E, A) <f (r, E, A) + tW s (E, A) g s (r, E, A) , (7) 

The superscripts v and s refer to volume and surface peaked terms respectively. The "COSH" 
form 4] is used in this paper, viz. 

{cosh[R(E,A)/a(E,A)]-l} 



{cosh [R(E, A) /a(E, A)] + cosh [r/a(E, A)] — 2}' 

= {cosh[R(E,A)/a(E,A)] - 1} {cosh [r / a(E, A)] - 1} 
{cosh [R(E, A)/a(E, A)] + cosh [r/a(E, A)] - 2} 2 

Similar forms have been taken for g v and g s . Thus the optical potential consists of scalar 
and vector terms each having real and imaginary parts. In the global analyses |4j the eight 
strength parameters varied with the proton center of mass energy E (in MeV) and with the 
atomic mass of the target A. The results can be represented by polynomials, 

4 3 

V V (E, A) = v + v mX m + ^ v n+4V n + v 8 xy + v 9 x 2 y + v 10 xy 2 , (9) 

m=l n=l 

where x = 1000/i? and y = Aj (A + 20). We have used the same form for all other potentials 
V S (E,A), W V (E,A), and W S (E, A). But the scalar and vector potentials, and their real 
and imaginary parts, all may have different geometry parameter values, the energy and mass 
dependence of which can be expressed by 



r o + TmXm + X] rn +^y n + r8X y + r ^ %2 y + r ^ x v 1 

771=1 71=1 
4 3 



R = A 1/z 

a = a + Y a mx m + Y a ^y n + a ?> x y + a $ x2 y + a w x y 2 - ( 10 ) 



777=1 77=1 
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There are many parameters in this specification which have been defined suitably by using a 
large data set in a search process. In the global DP calculations we use the recoil correction 
given by Cooper and Jennings 15]. 



C. Phase shifts, 5-matrices, and observables 

Irrespective of how the NA optical potential is specified, the objective of its use is to 
define the scattering (S')-matrix, or equivalently the (complex) phase shifts 5^(k). With 
superscripts (±) identifying j = I ± 1/2, these relate by 

S±(k) = e 2iS ?W = v ±(k)e 2m i s ^ where rjf(k) = \S±(k)\ = e ~ 2 ^W] . (11) 

With E oc k 2 , integral observables for neutrons, namely the total elastic, total reaction 
(absorption), and total cross sections then are given by 

_ 00 

a ^ E ) = ^J2{( l + 1 )\ s ^( k )- 1 \ 2 + l \ s r(k)-i\ 2 } , 



z=o 

_ oo 

^ = + [i-vm 2 ]+i[i- vr (k) 2 ]} , 

1=0 

and 

cttot(E) = a el (E) + a R (E) 
2?r °° 

= T2 E W + X ) I 1 - cos (2K [*+(*)])] r,r(k) cos (2» (*)])] } (12) 

respectively. Note that the Coulomb interaction means that only total reaction cross sections 
(formed by a modification of the above) can be measured from proton scattering. 

To evaluate angular observables, and specifically the differential cross sections and ana- 
lyzing powers, one need form the scattering amplitudes which, for nucleon elastic scattering 
are then 2x2 matrices in the nucleon spin space having the form, 

f(6) = A(6)+iB(6)tr-n, (13) 

where, in terms of the S-matrix elements 



1=0 

B ^ = -gk E - ^rc*)] plifi) . (M) 
i=i 

From these (complex) amplitudes, and with the quantization axis normal to the scattering 
plane, the (elastic scattering) differential cross section, analyzing power, and spin rotation 
are defined by 

<h \A(9)\ 2 + \B(e)\ 2 



dfl 

Ay{9) = 



23 \A*(d)B(d)} 
da/dQ 



6 



The g-folding model results have been obtained using the DWBA98 code [13| which 
evaluates scattering amplitudes using a helicity formalism Those scattering amplitudes 
A(6),B(6) relate to A(6),B(6) by a unitary transformation involving the rotation matrices 
{da(^)} under which the cross section and analyzing power are invariant, in that 

' h \A(6)f + \B(9)\ 2 



dtt 

Ay{9) = 

but the spin rotation becomes 
1 



A(9) + B{6) 



23 



A*{8)B{6) 



da /dVL 



Q(0) 



da /dVL 



A{6) - B{6) sin(0) 



A\6)B{e) cos(#) 



(16) 



(17) 



III. RESULTS AND DATA COMPARISONS 



We have chosen five nuclei for study They are 12 C, 16 0, 40 Ca, 90 Zr, and 208 Pb. Not only 
do they span a large range of target mass and ground state isospin, but also some have been 
described by mean fields derived from large space structure studies which suit comparison 
of scattering models. 

All microscopic model results we show have been evaluated using the DWBA98 pro- 
gram 0], input to which are the density dependent and complex effective NN interactions 
described earlier. Other input to DWBA98 are the ground state occupancies (OBDME) and 
the associated SP functions; both of which are defined from nuclear structure calculations. 
The structures we have used for each nucleus in the (/-folding calculations are 

• 12 C: For this, the lightest mass nucleus we consider, a no-core shell model calcula- 
tion [TtJ in the full (0 + 2)Huj space defined its ground state. The WBT interaction of 
Warburton and Brown [l8| was used. With exception of the 0^~ state, this shell model 
calculation gave a very good spectrum to 20 MeV excitation. 

• 16 O: For this nucleus, we have used the structure as determined by Haxton and John- 
son 0]. That was a (0 + 2 + 4)hw calculation made using a hybrid interaction based 
upon the Op-shell Cohen and Kurath |20J] potentials. The (positive parity) spectrum 
of 16 to about 10 MeV excitation is quite well reproduced. 



40 Ca: The ground state structure for this nucleus was taken from recent studies |7L l2lj 
made using the SHF(SKX) model of structure. That model lead to better fits to the 
200 MeV data (which we consider again herein) than either the simple Ohw shell model 
or SHF results made with other interactions. 



90 



Zr: For this nucleus, a shell model space defined in the 1»3, fs , lpi , Ogs orbits was 

■ 1 2 2 2 2 

used with the NIS interaction [22j. That interaction was optimized to describe nuclei 
in the Ni— >Sn mass region. 

208 Pb: This ground state of this nucleus also is described by a SHF calculation. In this 
case, using the SKM* interaction [?} , wave functions and OBDME (shell occupancies) 
result with which the 40, 65, and 200 MeV proton elastic differential cross sections as 
well as the electron scattering form factor have been extremely well reproduced |2l| . 
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As noted before |2|, with all details preset, the (/-folding approach is predictive. So only one 
calculation of an optical potential and of its use to give cross sections and spin observables is 
made. How good a reproduction of data results then is a measure of the formulation and/or 
of the specifications used. 

The global Dirac approach [4j is different in philosophy to the above. First it is phe- 
nomenological and the functional forms of parameter values were obtained by fitting a large 
body of data, including most if not all of the data considered herein. From those studies, two 
parameter sets, which we identify as the EDAI (^-dependent A-independent) and ED AD 
(independent A-dependent), specify the Dirac potentials we shall use. From the general 
form of the model given in Eq. (j7J), there are eight potentials which are functions of only 
energy and radius. The ED AD fits test the sensitivity of the computation to the input 
optical model potentials; though there are three different ED AD fits that give equally 
high-quality agreement with the data. 

In all figures to follow, results found using the (/-folding, the EDAI , and the EDAD3 
models will be displayed by the solid, long-dashed, and dot-dashed curves respectively. 



A. Cross sections and spin observables: angular variations. 

The differential cross sections for 65 MeV and 200 MeV protons elastically scattered 
from 12 C, 16 0, 40 Ca, 90 Zr, and 208 Pb are shown in Figs. and ^respectively. The calculated 
results are compared with data in nine of the ten cases. The data are displayed by different 
symbols: filled circles for 12 C, open squares for 16 0, open diamonds for 40 Ca, crosses for 
90 Zr, and open circles for 208 Pb. With all results and data scaled for clarity, from the top 
down in Figs, dare the results and data for 12 C (1), 16 (10" 1 ), 40 Ca (10~ 2 ), 90 Zr (lO^ 3 ), 
and 208 Pb (10~ 5 ). The scales used are given in brackets following the target symbol. 

Clearly both global DP results fit the data very well as the studies to determine the 
global set were designed to do. The (/-folding results, by comparison, do not do as well. But 
recall that each (/-folding result is a prediction built upon a given microscopic model of the 
target structure and an effective interaction determined from infinite matter solutions of the 
NN (/-matrices. As such, these results also are very good, predicting reasonable structures 
and essentially the correct magnitudes for all cases. The structures of the (/-folding results 
are sharper than those of the data suggesting that the absorptive character of the (/-folding 
optical potentials at 65 MeV may be a little weak. 

Such is far less the case with the 200 MeV results which are displayed in Fig. EJ In this 
figure the notation is as for Fig. ^ with now the scale for the 208 Pb results (lowest in the 
diagram) being 10~ 4 . Again both of the global DP cross sections give very good fits to the 
data while now the (/-folding predictions are good representations of that data though with 
some noticeable disparities at the larger scattering angles; where the magnitudes of the cross 
sections are less than 0.1 mb/sr and the momentum transfer values are in excess of ~ 2.5 
fm -1 . For such conditions, the first order theory upon which the (/-folding model is built 
may be too limited Q. 

The comparisons between data and calculated results for the spin observables are dis- 
played in Figs. EJ IH El and El In the first two results for analyzing powers A y are displayed, 
while in the latter, two spin polarizations Q are presented. The results for 65 MeV and 
200 MeV are shown on the left and right respectively in each figure. The available data are 
presented by the diverse symbols in each diagram while the calculated results are displayed 
by the curves. The notation is as specified previously. In Fig. El the analyzing powers from 
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FIG. 1: (Colour online ) T he differential cross sections from the scattering of 65 MeV protons 



from five nuclei. Data |23|, |2J, [25j are compared with the results of different theoretical model 



calculations. Details are given in the text. 



elastic scattering of 65 MeV and 200 MeV protons from 12 C, 16 0, and Ca are shown. Those 
from 90 Zr and 208 Pb are presented in Fig. HJ Once more the global DP calculations lead to 
very good fits to most of these data. 

Disparities between the EDAI and EDAD3 results, and between those results and data, 
are evident at large scattering angles for the 12 C and 16 cases. But, as noted before, those 
regions involve high momentum transfer and very small cross-section values; features that 
are hardest to fit within a global parameter scheme, let alone with a g-folding approach. 
Of course, if each target is treated independently then it should be possible to find better 
phenomenological fits to these higher momentum transfer data. As with the cross sections, 
the predictions obtained with the (/-folding approach, while not reproducing the data as 
well as the global DP, still are quite good in comparison with the data. In some cases, the 
g-folding predictions are very good fits. Overall then, the ingredients used in the (/-folding 
method can be considered quite realistic. 

The spin rotations Q resulting from the three calculations are compared in Figs. Eland El 
There are very little data available, but that for 208 Pb is shown. All three calculated results 
agree very well with that (limited) data set. For the other nuclei, the global DP results are 
quite similar, diverging from each similarly to the differences between the analyzing power 
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FIG. 2: (Colour online) The differential cross sections from the scattering of 200 MeV protons 
from five nuclei. Data are compared with the results of different theoretical model 

calculations. Details are given in the text. 

results. They also differ from the g-folding results in like fashion. As with the comparisons of 
analyzing powers, the spin rotations from all three calculations are in quite good agreement 
structurally to over 40° in the center of mass. 

B. Total reaction and total cross section results 

1. Proton total reaction cross sections 

In this subsection predictions of the total reaction cross sections for proton scattering 
are displayed. Results for energies to 800 MeV for the five nuclei considered are compared 
against data; the latter taken from many sources as listed in Table HI 

The results for scattering from 12 C are presented in Fig. [7| Total proton reaction cross 
sections are shown for a projectile energy range from 20 MeV to 900 MeV. Although ex- 
perimental data exist to much lower energies, any effective mean field prescription for the 
optical potential at such low energies is not appropriate. However, in a recent paper [52j], a 
multi-channel algebraic scattering (MCAS) theory has been developed with which the low 
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FIG. 3: (Colour online) Analyzing powers from 65 (left) and 200 (right) MeV proton elastic 
scattering from 12 C (top), 16 (middle), and 40 Ca (bottom). Details are given in the text. 

TABLE I: Sources of proton reaction cross-section data. 



Nucleus Proton references (in year order) 

[22], [30J, [31], [22], [33], [34J, [35], [36], [37J, [38], [39], [40], [41], [42], [43] 
16 [44], [4pJ, [45] 

40 Ca [34], [46], [38], [47], [45J, [43] 

90 Zr [48], [38], [49], [39] 

208 Pb [30], [32], [33], [35], [46], [50], [38], [39], [40], [51], [45], [43] 
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FIG. 4: (Colour online) Analyzing powers from 65 (left) and 200 (right) MeV proton elastic 
scattering from 90 Zr (top) and 208 Pb (bottom). Details are given in the text. 

energy range should be covered. This MCAS approach lends itself to specification of the 
dynamic polarization potential resulting from a summation over all diverse channel func- 
tions allowed in that theory. But the resultant optical potential then is extremely non-local 
and energy dependent; as it need be not only to produce a smooth energy varying elastic 
scattering cross section but also to give the sharp and broad resonances readily observed 
at low (< 20 MeV) energies (52 ]. The low level densities in most spectra to ~ 20 MeV is 
not a condition conducive to confidence in the (/-folding or any mean field optical model 
potentials. Nonetheless for 12 C, there are many data points at energies between 20 and 40 
MeV and a significant number thereafter to 500 MeV. Above that there are only a few. 

In Fig. Q the data to 300 MeV are well reproduced by the (/-folding model results. But 
for the higher energies, the ^-folding approach is inadequate; seriously under-predicting the 
data. That is the case for all five targets emphasizing the energy regime limit one need 
remember for this microscopic Schrodinger approach to scattering. On the other hand the 
reaction cross sections determined with both the EDAI and EDAD3 global potentials 
reflect the data very well falling within the uncertainties at most energies. Some data, 
notably at 61 MeV |33[ and at 77 MeV (35| are in disagreement with all three calculated 
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FIG. 5: (Colour online) Spin rotations from 65 (left) and 200 (right) MeV proton elastic scattering 
from 12 C (top), 16 (middle), and 40 Ca (bottom). Details are given in the text. 

results. However, they are exceptional poin ts in that they also disagree with empirical values 
taken at nearby energies. Menet et al. [39] argue that a much larger systematic error should 
have been used in the data processing of earlier experiments that lead to the exceptional 
points. Such also occur with data taken on other targets. 

Predictions for proton scattering from 16 are compared with data in Fig. |HJ There 
is no p- 16 total reaction cross section data for energies between 50 MeV and 200 MeV. 
Nonetheless the g-folding potential predictions to 50 MeV are in very good agreement with 
data for those energies. But as with the case of 12 C, the ^-folding approach under-predicted 
the four data points in the energy range 200 MeV and 600 MeV. These high energy data are 
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FIG. 6: (Colour online) Spin rotations from 65 (left) and 200 (right) MeV proton elastic scattering 
from 90 Zr (top) and 208 Pb (bottom). Details are given in the text. 

well replicated by the global DP potentials however, and by the EDAD3 model in particular. 
That model also gave results that compare well the data at lower energies while the EDAI 
result over-predicted them. 

Predictions for proton scattering from 40 Ca are compared with data in Fig. El There are 
relatively few data points in the energy range with but one (at 700 MeV) above 180 MeV. 
Again the predictions made using the ^-folding potentials are in very good agreement with 
the data for energies to 180 MeV. But the microscopic approach underestimates the datum 
at 700 MeV. Both the global DP (EDAI and ED AD?,) models replicate the data well at 
the energies to 65 MeV, overestimate the data at the energies near 100 MeV and 180 MeV, 
and reproduce very well the datum at 700 MeV. 

In Fig. El we present the data and our predictions of the total reaction cross sections for 
proton scattering from 90 Zr. There is very little data with which to compare, nevertheless 
the g-folding potential gives results in very good agreement with the five actual data points 
in the energy range 30 to 100 MeV. On the other hand both global DP model calculations 
slightly overestimate that data set. Again however, it is in the higher energy ranges that 
the ^-folding and Dirac model predictions diverge, as we have come to expect. 
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FIG. 7: (Colour online) Energy variation of ctr for proton scattering from 12 C. The curves are 
identified in the text while data references are as given in Table QJ 

In Fig. ^2 we compare the calculated total proton reaction cross sections from 208 Pb 
with data taken from the references listed in Table |1J For this nucleus, the data are quite 
scattered with some being quite disparate. In this case the ^-folding model calculations 
underestimate data at energies between 20 MeV and 40 MeV but agree well with the data 
taken at energies thereafter to 300 MeV. On the other hand, both EDAI and EDAD3 global 
DP models reproduce a select set of the data over the whole energy range quite well. As 
with 12 C, exceptional data points are found near 30 MeV, 61 MeV and 77 MeV. In this case, 
other data from 208 Pb taken at 60.8 MeV (3^ and 65.5 MeV |43[ give different results and 
in fact are values of reaction cross sections that are in accord with our model predictions. 

2. Energy variations of neutron total cross sections 

Accurate measurements of neutron total reaction cross sections are far more difficult to 
achieve than their proton counterparts. Indeed usually those cross sections are obtained by 
subtracting the elastic from the total scattering cross section. While both of those cross 
sections can be measured with some accuracy, the subtraction of two large numbers with 
attendant uncertainties compounds uncertainty in the result. We note, however, that a new 
Japanese technique, for measuring neutron total reaction cross sections utilizing in-beam and 
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FIG. 8: (Colour online) Energy variation of ctr for proton scattering from 16 0. The notation is as 
in the text. 

out-beam methods similar to those used in proton scattering, shows promise. Nevertheless, 
herein, we concentrate on analyses of total scattering cross-section data. 

The data we have chosen to analyze have been taken from a recent survey by Abfalterer et 
al. [53j] . That survey includes data measured at LANSCE that are supplementary and 
additional to those published earlier by Finlay et al. 54]. For comparison with that recent 
data compilation, we calculate neutron total scattering cross sections for the five nuclei of 
interest; 12 C, 16 0, 40 Ca, 90 Zr, and 208 Pb . 

In forming neutron optical potentials by the g-folding method, we have used the same 
structure models and effective NN interactions chosen to create the proton optical potentials 
whose results were presented and discussed above. Likewise global DP model potentials for 
neutrons were formed using the EDAI and EDAD3 parameter sets that have been used to 
get proton scattering results, and with the Coulomb potential set to zero, generated results 
displayed and discussed. No Coulomb energy shift has been made in those calculations. 
Thus all of the results shown in Fig. ^Jare predictions. For 12 C, neutron total cross sections 
are well predicted by all three model calculations. The EDAI and EDAD3 potentials of 
the global DP give very good results for energies above 20 MeV, while the ^-folding model 
calculations do so for energies between 50 MeV and 300 MeV. Outside of that range the 
g-folding results under-predict data. Below 20 MeV in this case, as well as for 16 which 
is considered next, resonance effects are very evident. Such are not encompassed in the 
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FIG. 9: (Colour online) Energy dependencies of ctr for proton scattering from 40 Ca. The notation 
is as in the text. 



mean field models we consider. For 16 0, the model results follow a similar pattern as for 
12 C. Global DP results give good fits to the data as does the ^-folding model but only for 
energies above ~ 70 MeV. The ED ADZ model predicts well data for the whole energy range 
while the EDAI model slightly over-predicts data between 20 MeV and 60 MeV. 

For the three heavier nuclei, 40 Ca, 90 Zr, and 208 Pb, the total cross sections down to 
10 MeV do not show significant resonance features but there is the onset of a Ramsauer- 
like effect which accentuates with target mass. None of the models predict such large scale 
oscillations adequately, though with improved structure details, as used in the calculations 
for scattering from 208 Pb, those oscillations are more noticeable when compared to results 
from use of a simple oscillator model of structure [Hij]. The results found for these three 
nuclei show a degradation in the quality of the predictions from all three model calculations. 
For 40 Ca, the comparisons are still similar to what is noticed for the lighter mass targets, 
though the g-folding model predictions above 300 MeV are poorer. For 90 Zr and 208 Pb all 
three model results clearly depart from what is observed. However, the g-folding results 
trend quite well as average results upon which a Ramsauer effect may be superimposed |56| . 
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FIG. 10: (Colour online) Energy dependencies of an for proton scattering from 90 Zr. Basic notation 
is as for Fig. [7| 

IV. CONCLUSIONS 

A microscopic coordinate space Schrodinger model of nucleon-nucleus interactions and 
global DP model for the same systems, have been used to predict proton and neutron 
scattering cross sections; both angular and integral. Cross-section data and spin observables 
have been well described for energies at which the diverse models are valid. That includes 
a substantial range of energies in common for both approaches. 

Our first study was of angular dependent observables; differential cross sections and spin 
observables (analyzing powers and spin rotations) for elastic proton scattering at 65 MeV 
and at 200 MeV. Both methods of analysis gave good reproduction of measured data. The 
two models, with the global DP potentials defined from the parameter specifications in the 
recent study of their energy and mass dependences were used to predict proton total 
reaction and neutron total cross sections in scattering from five nuclei, 12 C, 16 0, 40 Ca, 90 Zr, 
and 208 Pb. Both models gave good predictions of that data (in the energy regimes for which 
the models are appropriate) with the largest discrepancies reflecting the Ramsauer effect in 
neutron total cross sections. 

Fundamentally, our Dirac and Schrodinger approaches are different. The microscopic 
Schrodinger model incorporates all the dominant medium modifications in the optical po- 
tential, without significant approximation, by using a realistic ground state density to give 
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FIG. 11: (Colour online) Energy dependencies of ctr for proton scattering from 208 Pb. The notation 
is as in the text. 

a reasonable specification of all terms in the optical potential. The global DP approach 
provides a natural specification of such terms in local equivalent form. Nonetheless, the 
comparisons obtained in this work are exemplary of the dilemma in judging the relative 
merits of relativistic vs. non-relativistic approaches for the analyses of intermediate energy 
nucleon-nucleus scattering. That the two models give results that agree so well for the dif- 
ferential cross sections as well as for both spin observables and for both projectiles gives 
confidence in the reality of calculated results at the two energies we considered. As has been 
speculated 0] however, the answer throughout the energy regime may lie in QCD-based 
models of nuclear scattering systems. Indeed the work of Cohen et al. [57] shows evidence 
that the large scalar and vector fields of Dirac Phenomenology may be related to quark 
degrees of freedom in the nucleon. However, such QCD-based theoretical models and con- 
comitant experiments will require simultaneous treatment of proton and neutron scattering 
in order to be complete. Sadly, the latter's experimental database currently is lacking. 
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FIG. 12: (Colour online) Total cross sections for neutrons scattered from 12 C, 16 0, 40 Ca, 90 Zr, and 
208 Pb. The data are those of Abfalterer el al. [H^]; for 12 C the data correspond to scattering from 
natural carbon. 
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